UH-511-1123-08 
hep-ph/yymmnnn 



0{a 2 s ) corrections to semileptonic decay b 

Kirill Mclnikov 



Department of Physics and Astronomy, 
University of Hawaii, 
Honolulu, HI, D96822 



00 

o 
o 

(N 



Oh' 

43 . 



in . 
Ch ■ 

q ; 

co ■ 

o : 

00 . 

o ■ 



X 



We compute the next-to-next-to-leading order (NNLO) QCD corrections to b —* clv\ decay rate 
at a fully differential level. Arbitrary cuts on kinematic variables of the decay products can be 
imposed. Our computation can be used to study the NNLO QCD corrections to the total decay 
rate as well as to the lepton energy, hadronic invariant mass and hadronic energy moments and to 
incorporate those corrections into global fits of inclusive semileptonic B-decays. 
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I. INTRODUCTION 



Inclusive semileptonic decays of _B-mesons into 
charmed final states are benchmark processes at B- 
factories. Because of relatively large rates and clean ex- 
perimental signatures, these decays can be studied with 
great precision. On the other hand, theoretical descrip- 
tion of semileptonic B decays is robust thanks to the 
Operator Product Expansion (OPE) in inverse powers 
of the 6-quark mass nib- The application of the OPE 
to semileptonic decays of _B-mesons leads to the con- 
clusion that both the total decay rate and various kine- 
matic distributions can be described by power series in 
Aqcd /TUb [H- For infinitely heavy 6-quark, the decay 
rate B — > X c lvi coincides with the rate computed at the 
quark level. For realistic values of bottom and charm 
masses, a few non-perturbative matrix elements that en- 
ter at order (Aqcd / Vra&) n , n = 2,3 are accounted for in 
existing theoretical predictions. 

In recent years, many measurements of moments of 
charged lepton energy and hadronic invariant mass in 
B — > X c Wi decays have been performed by BABAR, 
BELLE, CLEO, CDF and DELPHI i, i, i, . Com- 
parison of these experimental results with theoretical pre- 
dictions for corresponding observables leads to the de- 
termination of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element | V^t, | , bottom and charm quark masses 
and a number of non-perturbative parameters such as (j% 
and (1q 0, 0] ■ Typical precision claimed in these analyses 
is about one percent for |V^b| an d nib and a few percent 
for m c and non-perturbative matrix elements 0, 0] . 

To achieve such precision, advances in theoretical un- 
derstanding of semileptonic -B-decays were necessary, in- 
cluding subtle interplay between perturbative and non- 
perturbative physics and significant developments in 
technology of multi-loop computations. While one-loop 



corrections to both b — > clvi total decay rate and a 
number of important differential distributions are known 
since long [10|], it is interesting to remark that phc- 
nomenologically relevant triple differential distribution 
in charged lepton energy, leptonic invariant mass and 
hadronic invariant mass was computed through 0(a s ) 
only a few years ago [ll|, [l2j]- This fact illustrates the 
complexity of perturbative calculations, when massive 
particles are involved, at a fully differential level. 

Given the precision of available experimental measure- 
ments, good understanding of non-perturbative effects 
and a fairly large value of the strong coupling constant 
a s (m,b) « 0.24, it is expected that 0{a 2 s ) corrections to 
b — ► X c Wi decays are required for a consistent theoretical 
description. However, as was realized long ago, technical 
complexity of such an endeavor is daunting. To simplify 
the problem, the 0(p? s ) corrections to b — > X c lvi were 
computed in three specific kinematic points [HI, [HI, [l(| . 
These results were used in Ref . [l3[ to estimate the NNLO 
QCD corrections to T(b — ► X c li>i). Unfortunately, such 
a description is necessarily limited in its scope even for 
the total rate and a generalization of such an approach 
to more differential quantities, such as lepton energy and 
hadronic invariant mass moments, is clearly out of ques- 
tion. 

On the other hand, a subset of the NNLO QCD cor- 
rections, the BLM corrections [T3|, received significant 
attention recently. The BLM corrections are associated 
with the running of the strong coupling constant; they 
are potentially important since the QCD /3-function is 
large. For B decays, however, the BLM corrections are 
known to be modest if proper definition of quark masses 
is adopted and judicious choice of the renormalization 
scale in the strong coupling constant is made. The BLM 
effects are the easiest NNLO effects to calculate since 
they can be obtained from a one-loop computation if 
the latter is performed with a non- vanishing gluon mass 
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[1 81 ] . For this reason, in the past, the BLM corrections 
to b — > clvi were calculated for the total rate and various 
kinematic moments [H[2Qj|. However, the NNLO QCD 
corrections beyond the BLM approximation, for which 
genuine two-loop computations are required, remained 
missing. 

Calculation of these two-loop corrections became pos- 
sible recently thanks to developments in numerical ap- 
proaches to multi-loop computations [21| . These numer- 
ical methods benefit from the absence of mass hierarchy 
in the problem which is the case for b — > c decays, since 
masses of bottom and charm quarks are close. The possi- 
bility to use the approach of Ref . [2l[ to describe decays of 
charged particles was recently pointed out in [22j where 
electron energy spectrum in muon decay was computed 
through second order in the perturbative expansion in 
QED. 

The goal of this Letter is to present the computation 
of 0{a 2 a ) corrections to b — > XJPi decay rate at a fully 
differential level. Our results can be used to calculate ar- 
bitrary observables related to inclusive 6 -t c transition 
through NNLO in QCD. For example, second order QCD 
corrections to such popular observables as lepton energy, 
hadronic invariant mass and hadronic energy moments 
can be studied in dependence of the cut on charged lep- 
ton energy. Inclusion of the results of our computation 
into global fits, should lead to a reduction of the theoreti- 
cal uncertainty in the determination of | V^b | , the bottom 
and charm quark masses and the non-perturbative pa- 
rameters that contribute to the decay rate. 



II. COMPUTATION 

In this Section, we set up our notation and briefly de- 
scribe technical aspects of the computation. A detailed 
description of the method can be found in [U [22j • 

Consider the decay b — > XjDi where the final state 
lepton is massless. The differential decay rate can be 
written as 



dr 



G F \V ch \ 2 ml 
192vr 3 



(dfb + a dFi + a 2 dF 2 



(1) 



where Gf is the Fermi constant, nib is the 6-quark pole 
mass, a — a s /ir and a s is the MS strong coupling con- 
stant defined in the theory with five active flavors and 
renormalized at the scale mj. For numerical computa- 
tions, we use mt = 4.6 GeV and m c = 1.15 GeV. While 
these numerical values for the quark masses can not be 
justified in the pole scheme, our choice is motivated by 
an eventual necessity to transform the pole scheme com- 
putation to a more suitable scheme. The values of the 
quark masses that we employ in this Letter correspond 
to the central values of m^ c in the "kinetic scheme" [23j j . 
derived in recent fits to inclusive semileptonic 73-decays 
0,|. 

To calculate the functions <1Fq_2, we have to account 
for different processes. At leading order, (IFq is computed 



by squaring the matrix element of the process b — > cWi 
and summing or averaging over spins and colors, as ap- 
propriate. At next-to-leading order, d7<i receives contri- 
butions from virtual 0(a s ) corrections to b — > clvi and 
from the real-emission process b — ► clDi + g. To compute 
d7<2, we require two-loop 0{a 2 s ) corrections to b — > clT>i, 
one-loop 0(a s ) corrections to b — > cli>i+g and the double 
real-emission corrections b — > clui+X, where X refers to 
two gluons or a quark-antiquark pair or a ghost-antighost 
pair. We will refer to these corrections as double-virtual, 
real-virtual and double-real, respectively. In addition, 
we have to account for a variety of renormalization con- 
stants, when computing higher order corrections. We do 
not include the process b — > cWi + cc in our calculation 
since the energy release in this process is so small that it 
can not be treated perturbatively. 

To calculate the NNLO QCD corrections, the method 
for multiloop computations developed in [2l|, [12] is em- 
ployed; in those references a detailed discussion of many 
technical issues relevant for the current computation can 
be found. One technical aspect that we improve upon rel- 
ative to Refs. [U, is how virtual corrections to single 
gluon emission process b — > cWi + g are treated. In Refs. 
pH |22| these corrections were dealt with by an analytic 
reduction to master integrals followed by a numerical 
evaluation of those. This method, however, becomes im- 
practical quite rapidly, once the number of external par- 
ticles or the number of massive particles in the problem 
increases. In principle, the real- virtual corrections can be 
computed numerically, but for heavy-to-light decays this 
is complicated because some Feynman diagrams develop 
imaginary parts. To handle these imaginary parts, we 
proceed as follows. For all Feynman diagrams that con- 
tribute to real- virtual corrections, it turns out possible 
to identify a Feynman parameter that enters the denom- 
inator of the integrand linearly. Let us call this Feynman 
parameter x\. Then, a typical integral that has to be 
computed reads 



7(0,1) = 



(-a + 611 +i0) 1+e ' 



(2) 



Here n> — l,b > a > and both a and b depend on other 
Feynman parameters and the kinematic variables. The 
two arguments of the function 7 refer to lower and upper 
limits for x\ integration. To calculate 7(0,1), we note 
that by extending upper integration boundary to infinity, 
a solvable integral for arbitrary a, b and n is obtained. On 
the other hand, since 



1(0,1) =7(0, 00) -1(1, 00), 



(3) 



and because the denominator of the integrand in 7(1, 00) 
is sign-definite, 7(1, 00) can be computed numerically in 
a straightforward way. It turns out that, up to minor 
modifications, this trick can be used to avoid dealing with 
the imaginary parts for all Feynman diagrams that con- 
tribute to one-loop corrections to single gluon emission 
process in b — > c decays. 
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Because couplings of quarks and leptons to the charged 
current are chiral, proper treatment of the Dirac matrix 
75 in d = 4 — 2e dimensions is important. While this 
problem can be avoided in the computation of the total 
decay rate, for more differential quantities it becomes an 
issue. We use the approach of Ref. [24[ where a consistent 
framework for extending the axial vector current to d- 
dimensions is given. 

Our computation can be checked in a number of ways. 
First, the double- virtual, real- virtual and double-real cor- 
rections are divergent when taken separately but these di- 
vergences must cancel in physical observables. We have 
checked these cancellations for a variety of observables, 
from the inclusive rate to various moments with cuts on 
both charged lepton energy and the hadronic invariant 
mass. Second, in the limit m c — » mj, the NNLO QCD 
corrections to the decay rate b — > clvi are described by 
the so-called zero-recoil form factors computed through 
0{o? s ) long time ago [l6[. We have checked that in the 
limit m c — > mj, our computation reproduces the zero- 
recoil form factors. Third, we can use published results 
for the BLM corrections to the total rate and charged 
lepton energy, hadronic invariant mass and hadronic en- 
ergy moments [20] to check parts of our computation re- 
lated to massless quark contributions to gluon vacuum 
polarization. Finally, considering the limit m c <C m,b, we 
reproduce the NNLO QCD corrections to b — ► uWi decay 
rate reported in Ref. [25| . 



lepton or hadronic system in the &-quark rest frame and 

The lepton energy moments introduced in Eq.(0]) can 
be written as 

L n = 4°) + aLM + a 2 (/?o4 2 ' BLM) + 4 2) ) + (7) 

where ellipses stands for higher order terms in the pertur- 
bative expansion in QCD. Similar decomposition can be 
performed for the hadronic energy moments H n . In ad- 
dition, we use Pq = 11 — 2/3Nf and define the non-BLM 

(2) (2) 

corrections L\ ' , H„ as the difference of the complete 
0{a 2 s ) correction and the BLM correction computed with 
N f = 3. 



n 


Scut, GeV 


r(0] 


rW 


r (2, BLM) 


r(2) 








1 


-1.77759 


-1.9170 


3.40 


1 





0.307202 


-0.55126 


-0.6179 


1.11 


2 





0.10299 


-0.1877 


-0.2175 


0.394 





1 


0.81483 


-1.4394 


-1.5999 


2.63 


1 


1 


0.27763 


-0.49755 


-0.5667 


1.00 


2 


1 


0.09793 


-0.17846 


-0.20875 


0.382 



TABLE I: Lepton energy moments. 



III. RESULTS 

We are now in position to discuss the results of our 
computation. We consider a number of observables, 
mostly for illustration purposes. We present the results 
in the pole mass scheme and use the strong coupling con- 
stant renormalized at the scale m&. While the pole mass 
scheme is known to be an unfortunate choice inasmuch 
as the convergence of the perturbative expansion is con- 
cerned, we decided to present our results in this way for 
clarity. However, we emphasize that the impact of the 
NNLO QCD corrections, computed in this paper, on the 
determination of \V c h\, heavy quark masses and the non- 
perturbative parameters, including kinetic and chromo- 
magnetic heavy quark operators, can only be assessed 
once the pole mass scheme is abandoned in favor of a 
more suitable quark mass definition and the NNLO QCD 
corrections are included into the fit. 

To present the results, we follow Ref. [2(| and define 

T (T? x ((Ei/m b ) n 6JA - gout) dT) 
L n (E cut ) = — , (4) 

„ ,j? x i(E h /m b r 9{E l - E cut ) dT) 
H n (E cut ) = — , (5) 



In Tables HHTl the results for lepton energy and hadronic 
energy moments with and without a cut on the lep- 
ton energy are displayed. The numerical accuracy of 
L^\H^ and Ll 2 ' BLM) , i# BLM) is about 0.1 - 0.2% 
whereas the numerical accuracy of , is about 
1 — 3%. It is possible to improve on the accuracy but 
this requires somewhat large CPU time. Nevertheless, 
for all practical applications the achieved numerical ac- 
curacy is sufficient. 

There are a few interesting observations that follow 
from Tables IIIIII Quite generally, the non-BLM correc- 
tions and the BLM corrections have opposite signs; given 
their relative magnitude and the value of /?o, it is easy 
to see that the 0{a 2 s ) corrections are about twenty per- 
cent smaller than what the BLM-based estimates suggest. 
The relative magnitude of the non-BLM and BLM cor- 
rections is largely independent of n and of whether the 
lepton energy cut is applied. 



It 


E cut , GeV 


H (0) 

1 1 11 


1 1 n 


rr(2,BLM) 
' ' n 


H (*) 
11 n 


1 


1 


0.334 


-0.57728 


-0.6118 


1.02 


2 


1 


0.14111 


-0.23456 


-0.2343 


0.362 



TABLE II: Hadronic energy moments. 



where (...) denotes average over the phase-space of all 
final state particles, Ei : h is the energy of the charged 



First row in Table U provides the NNLO QCD cor- 
rections to the total decay rate b — > cWi in the pole 
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mass scheme. Such corrections were estimated earlier 
in Ref . [l3T] . Note that in Ref.[l3T] the numerical results 
are given for the ratio of quark masses m c /mb = 0.3 
and also the BLM corrections are defined with Nf = 4, 
rather than Nf = 3. Calculating the non-BLM correc- 
tions for the set of parameters employed in [l3[, we find 
L k 1.73 which is to be compared with the estimate 

w 0.9(3), reported in [HG3]. 
The results of Ref. [l3[ were used in Ref. [26| to estimate 
the impact of the QCD corrections on T(B — > X c lv{). In 
Ref. [261] the perturbative corrections to b — ■> clvi decay 
rate are described by a factor A pert , defined as 

T{b^X c Wi) = A pCTt (r) (dTo), (8) 

where r = m c jm] 3 . A pcrt depends on the adopted 
scheme for the quark masses. In the kinetic mass scheme, 
AP° rt (0.25) = 0.908 is quoted. To arrive at this re- 
suit, Ref. [26j uses L — 1.4 which is about a factor 2.5 
smaller than the corresponding entry in Table HI Cor- 
recting for this discrepancy, we derive 

A port (0.25) = 0.919. (9) 

We believe that this value for the perturbative renormal- 
ization factor in the kinetic scheme for m c /mi, = 0.25 
should be employed in global fits of semileptonic B- 
decays. 

Further analysis of entries in Table [T] suggests that the 
QCD corrections in general and the non-BLM correc- 
tions in particular mostly affect the overall normalization 
rather than shapes of kinematic distributions. This fol- 
lows from the approximate independence of Ln' 2 ^/Ln 
of n and also of whether or not the cut on the lcpton 
energy is imposed. It is therefore possible to speculate 
that the non-BLM corrections computed in this Letter 
will mostly affect the extraction of | V c h | whereas their in- 
fluence on, e.g., the 6-quark mass determination will be 
minor. Concerning the | V^t, | , the increase of the pertur- 
bative renormalization factor ^4 por t by 10 x 10~ 3 implies 
the change in the value of |Kb|j extracted in Ref.Q], by 
about —0.25 x 10~ 3 . On the other hand, since non-BLM 
corrections were not included in a fit of Ref. [?J , the shift 
in the value of \V C \,\ derived in that reference will likely 
be larger ~ —0.5 x 10^ 3 . Although expected shifts in 
central values of |Vct>| are n °t large , we stress that they 
are comparable to uncertainties in |T4b|, derived in the 
global fits 0,1. 



IV. CONCLUSIONS 

In this Letter, the computation of the NNLO QCD 
corrections to the fully differential b — > clvi decay rate 
is reported. The differential nature of the computation 
makes it possible to apply arbitrary cuts on the kine- 
matic variables of final state particles. This result allows 
to extend the existing determinations of the CKM ma- 
trix element |T4b|, the bottom and charm quark masses 
and the non-perturbative parameters fj% and [1q from 
global fits to semileptonic decays of B-mesons, by in- 
cluding the NNLO QCD corrections exactly. We note 
that for a consistent high-precision analysis of semilep- 
tonic S-decays, also 0(a s ) corrections to Wilson coef- 
ficients of non-perturbative kinetic and chromomagnetic 
operators are required. Such a correction is available for 
the kinetic operator [27J but is still missing for the chro- 
momagnetic. 

We presented a few results for charged lepton energy 
moments and hadronic energy moments with and with- 
out a cut on the lepton energy in the pole mass scheme. 
These results suggest that the magnitude of the non-BLM 
corrections does not depend strongly on the kinemat- 
ics; the non-BLM corrections are approximately 2% for 
all the moments considered. We therefore expect that 
the non-BLM NNLO QCD corrections will mostly affect 
the determination of \V c b\ decreasing its central value by 
about one percent whereas their impact on the quark 
masses and the non-perturbative parameters will proba- 
bly be quite mild. 

As a final remark, we note that it would be interesting 
to extend this calculation in two ways. First, one may 
consider semileptonic decays of i?-mesons into massive 
leptons. Such an extension, relevant for the description 
of B — > X c + t + D T decay is straightforward. Second, it 
is interesting to extend the current calculations to allow 
for a massless quark in the final state. This is a difficult 
problem but it is highly relevant for the determination of 
the CKM matrix element |Kib| from semileptonic b — > u 
transitions. 
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